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adjustable constant 
computational domain 

integration domain size used for averaging (eq. (29b)) 
unit vectors in x- and y-directions, respectively 

reference height of background shear flow (reference length scale) 

size of computational domain in streamwise and vertical directions, 
respectively (fig. 2) 

of the order of 

Reynolds number of background flow, U^L/v 

distance from kth vortical center 

time 

2 

initial age of vortical spot, 5^(0 )/4v 


ups tr earn ( s panwi se ) f low 
reference velocity scale, Uq(°°) 

velocity components in x- and y-directions, respectively 
velocity vector 

circumferential velocity around kth vortical spot 
kth vortical center 
X^(t),Y^(t) velocity of kth vortical spot 


X R - x L horizontal spread of right and left vortex of pair 
x,y,z spanwise, vertical, and flight directions, respectively 
Y r vertical position of right vortex of pair 

T strength of vortical spot 

$ effective core size of vortical spot (small length scale) 

e small parameter, <$/L << 1 

C(x,y, t) vorticity distribution 
v kinematic viscosity 
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2 2 1 / 2 

distance from origin of coordinates, (x + y ) ' 

\J>(x,y,t) stream function 

ai(x,y,t) background vorticity distribution 

At time step 

Ax, Ay, AL grid sizes 
2 

V Laplacian operator 

Subscripts: 

k kth vortical spot 

L left 

min minimum 

R right 

x - 8/8x 

y = a/9y 

0 reference initial quantity at t = 0 

Notation: 

(~) bar quantity denotes terms of small parameter 

(^) tilde quantity denotes variation of background flow quantity 

( ) 1 derivative of quantity 

< > average of flow quantity (eq. (29b)) 
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SUMMARY 


The drift of trailing vortices in a crosswind near the ground is simulated by an 
unsteady, two-dimensional, rotational flow field with a concentration of large 
vorticity in "vortical spots" (having a finite but small effective size and finite 
total strength) . The problem is analyzed by a combination of the method of matched 
asymptotic analyses for the decay of the vortical spots and the Euler solution for 
the unsteady rotational flow. Using the method of averaging, a special numerical 
method is developed in which the grid size and time step depend only on the length 
and velocity scales of the background flow and are independent of the effective core 
size of a vortical spot. The core size can be much smaller than the grid size, 
whereas the peak velocity in the core is inversely proportional to the spot size. 
Numerical results are presented to demonstrate the strong interaction between the 
trajectories of the vortical spots and the change of the vorticity distribution in 
the background flow field. 


1 INTRODUCTION 

It has been well-recognized that there is a potential hazard whenever a small 
aircraft encounters the vortex wake of a larger aircraft, particularly during the 
takeoff and landing phase of the aircraft. An understanding of the motion and struc- 
ture of the vortex wake, including the ground effect and its minimization, is impor- 
tant for flight safety and also for efficient use of the airport (refs. 1 and 2). 

The variation of the vortex wake strength along the trailing edge of the lifting 
wing depends on the characteristics of the wing. Complete prediction of the flow 
field including the near wake is very difficult because of the disparate length 
scales associated with the generation, interaction, and eventual decay of the 
vortices. The problem of a steady far-wake vortex can be simplified by reducing the 
problem to an equivalent unsteady two-dimensional problem in a plane normal to the 
flight direction (ref. 3). This assumption ignores the streamwise curvature of the 
trailing vortex filaments, their initiation at the trailing edge, and the variation 
of the velocity parallel to the downstream direction (z) . Mathematically, it is 
assumed that 3/3z « 3/3x and 3/3z « 3/ 3y, where x and y are the s panwise and 
vertical directions, respectively. In the xy-plane at a station z, the trailing 
vortices are represented by "vortical spots" (having a small effective size inside 
which there is a strong vorticity distribution with finite total strength). 

This simulation is employed to study the drift and decay of far-field trailing 
vortices (vortical spots) in a crosswind (a spanwise shear flow) near the ground. 

Note that the Reynolds number R L of the background flow can be much larger than 1 , 
i.e.. 


R 


L 


V 


>> 1 


( 1 ) 



where U w , L, and v are the reference velocity, the height of the background shear 
flow, and the kinematic viscosity*, respectively. Since the size of a vortical spot 
is much smaller than L and its total strength T is finite, i.e., of the order 
of the vorticity and the velocity gradients in a spot can be very large. Thus, 

the diffusion in the core of a vortical spot is important; but away from it, the 
viscous term is of the order of 1/R L and, hence, is negligible. 

It is very inefficient to study the flow field of vortical spots submerged in a 
rotational flow by numerical solution of the Navier-Stokes equations for the entire 
flow field since the viscous terms are important only in small spots of high vor- 
ticity concentration and the grid size would have to be smaller than the size of the 
spots. A multiple-scale analysis is introduced in order to isolate the viscous decay 
of vorticity in the spots as a solution of the "small” scale variables, whereas the 
flow field in the "normal” scale obeys the Euler equations. The decay of the 
vortical core in each spot can then be described by a matched asymptotic solution 
(ref. 4), and the motion of the vortical spots is coupled with the variation of the 
vorticity distribution in the background nonuniform shear flow, requiring a numerical 
solution of the unsteady Euler equations. 

It is noted here that if the background flow is a constant shear flow (constant 
vorticity) or a potential flow (zero vorticity), the movements of a vortical spot 
cannot induce any variation in the background vorticity distribution and will not 
change the background flow. Therefore, the flow field can be represented as the 
superposition of the steady background flow and that induced by the moving vortical 
spots defined by the asymptotic analysis. 

Section 2 presents a description of the governing equations and the scaling laws 
for the mathematical simulation of the interaction of the vortex wake with a cross- 
wind. The asymptotic solutions for the decaying vortical spots are described. In 
this simulation, the background flow (without the vortical spots) fulfills the non- 
slip condition along the ground. The unsteady boundary layer along the wall (or 
ground) induced by a vortical spot has not been accounted for in the Euler solution. 
Therefore, the distance from a vortical spot to the wall should remain much larger 
than the thickness of the boundary layer induced by the vortical spots. 

In section 3, a special numerical method that implements the concept of a two- 
length scale analysis is described for the solution of the Euler equation. Conse- 
quently, the step size and the time step will be independent of the effective size of 
the vortical spots (which can be much smaller) . Note that the peak velocity in the 
core, being inversely proportional to the core size, can be much larger than U tt . 

In section 4, numerical results are presented to demonstrate the interactions 
between the vortical spot and the background nonuniform shear flow. The distinct 
features of the trajectories of vortical spots with varying strengths and different 
initial positions are illustrated. 


2 MATHEMATICAL FORMULATION 

2.1 The Physical Problem (Upstream and Initial Conditions) 

The physical simulation for the motion and decay of strong two-dimensional 
vortical spots interacting with a background rotational flow is illustrated in 
figure 1. An upstream (spanwise) flow U^(y) is specified, i.e.. 
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(2a) 


V(x,y,0) = u 0 (y) i = 0(1) 


as x + -«> for y > 0 with U q(°) = The ground is represented by y « 0. The 
background rotational flow has the reference length scale L, which is equal to the 
effective height of the shear flow, and the reference velocity scale U^, which is 
U-(°°). Both scales, L and U , have been set equal to 1, It is assumed that 

0 oo 

Uq(y) + 1 exponentially in y as y > <», e.g.. 


U Q (y) = 1 - e 


(2b) 


The corresponding initial background vorticity distribution is then 

w(x,y,0) * -U^(y) * 0(1) (2c) 

To simulate the trailing vortices, a highly concentrated vortical spot, say the kth 
spot centered at (X^(0), Y^(0)), is defined by the following characteristics: an 

initially small effective core size (^(°) << 1 ) and a strong vorticity distribution 
( = O(S^) >> 1 ) so that the total strength is finite ( : r - 0(1) = OdJ^L) ) . 

Because of these characteristics of vortical spots, it is assumed that the 
initial vorticity distribution £(x,y,0) can be split into two parts and that 


S(x,y,0) = u>(x,y,0) + e~ 2 £(x,y,0) 


(3) 


where 


e = 8 /L « 1 


(4) 


and 6 ^ is the typical core size* The second term e £(x,y,0) in equation (3) 

Mri'th C = 0(1) represents concentrated vorticity distributions near the vortical 
spot. This vorticity distribution is of compact support or decays exponentially in 
r^, where r^ is the distance from the kth center in terms of the small length scale 
5. One can then write £ as a function of the stretched variables x and y with 



for each k. For simplicity, the case considered here is that of vortical spots 
taving similar core structures (such as a Lamb vortex (ref. 5)), i.e.. 
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( 6 ) 


N 

-2- V -2 - - - 

£ C = ® ( x f y , 0 ) 

k=1 


N 


-z 

V=1 


r * -{yu/yon} 2 

2 e 

ir 6^(0) 


where 

? k = f[x - X k (0)] 2 + (y - Y k (0)] 2 } 1/2 /e 


From the initial vorticity distribution, it is noted that (1) the initial vortex 
positions (X^(0), Y^CO)) are assigned, (2) decays exponentially in r^, and 
(3) the strength T. is assumed to be of the order U L, i.e*, 

K 00 


oo 

-// 


e 2 l. dx dy 


00 

- If 


\ a \ a5 k * ° <a - L) 


(7) 


This analysis can take care of nonsimilar initial profiles, i.e., not of the 
type specified by equation (6), as long as £ decays exponentially in r and 
fulfills equation (7). 


2.2 Governing Equations 

For an unsteady viscous flow field, the vorticity distribution £(x,y,t) and 
the stream function i|/(x,y,t) are governed by the two equations 


H 

at 


+ -fe (ui) + ■ vv2t 


( 8 ) 


and 


* = 


-C 


(9) 


where 


V 2 


is the Laplacian operator. 


The vorticity is defined by 


_ j)v jto 
5 " 3x ' 3y 


( 10 ) 
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and the velocity components are related to the stream function by 


3i|> 3i p 

3 y 3 x 


( 11 ) 


The boundary conditions along the ground plane are 


u(x,0,t) = 0 


( 12 ) 


and 


v(x,0,t) = 0 


( 13 ) 


In general, to construct the Navier-Stokes solution of equations (8) and ( 9 ), one has 
to use grid sizes Ax and Ay << 6 and, hence, a very small At. This process is 
very inefficient since the viscous term vV 2 5 is important only in the neighborhood 
of a vortical spot (of the order 6); and away from it, the viscous term is not 
important. 

In section 3 , a special method is described that takes into account the special 
features of the vorticity distribution. It should be recognized that the flow field 
also has a small length scale 6, i.e., £L. The small parameter s is of the order 

of i,e *' 



so that the viscous terms remain important in the small scale, i.e., in the neighbor- 
hood of a vortical spot. On the other hand, in the normal scale, the flow field 
obeys the Euler equations. 


3 METHOD OF SOLUTION 

In subsection 3 . 1 , the problem is set up in two length scales, the normal scale 
L and the small scale 6 = eL. Analytical solutions describing the "small" scale 
flow field (i.e., the viscous decay of the strong vorticity in each vortical spot) 
will be identified as that of the similarity solution from the asymptotic analysis 
(refs. 6 and 7 ). The velocity (X^ft), Y^.(t)) of the kth vortical spot is coupled 
with the inviscid solutions describing the "normal” scale flow field (i.e., the back- 
ground flow). In subsection 3 . 2 , the details will be described about the method for 
constructing the numerical solution of the variation of vorticity in the background 
shear flow. A particular feature of the numerical method is that the grid size 
Ax = Ay = AL can be selected independent of the core size d . The time step 
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depends only on the grid size and on the length and velocity scales of the background 
flow, i • e . , 


At = f ( AL,U ) 
00 


(15) 


3.1 The Two-Length-Scale Problem 

The special form of the initial data suggests the following decomposition of 
vorticity: 

C(x,y, t) = eo (y) + C(x,y,t,e) + e 2 c(x,y,t,e) (16) 

Here, x and y denote the regular spatial variables and x ~ x/e and y = y/e 
denote the stretched variables. The first term ^w^(y) in equation (16) is the 
initial background vortiqity. The second term £(x,y,t,e) is the variation of back- 
ground vorticity induced by the vortical spots. The last term is composed of the 
large vorticity 0(e 2 ) near the vortical spot, i.e. , 

e~ 2 C = e ~ 2 (17) 

k 


Here, XL decays exponentially in 
x 


= {[x - X^t)] 2 + [y - Y k (t)J 2 } 1/2 /e 


Similarly, the stream function can be written as 


ij;(x,y,t) = ^q(y) + iKx,y,t, e) + ij/(x,y,t,e) (18) 


The velocity components can also be decomposed into 

V = [u Q (y) + u + e u]i + (v + e v)j (19) 
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and can be related to 


and ij> by 




and 


u 


ay 



( 20 ) 


( 21 ) 


The velocity fields are then described by the stream functions, which in turn are 
related to the vorticity functions by 


r 

= I o (y* ) dy' 
Jq 


( 22 ) 


V 2 4p = -l (23) 

and 

= -1 (24) 


The boundary condition on the ground (y - 0) for the velocity field is v(x,0,t) = 0. 
The nonslip condition shall be taken care of by the addition of a thin boundary layer 
induced by the vortical spots. By following the method of multiple scales, the de- 
composition equations (eqs._(16) to_(19)) are substituted into equations (8) and (9) 
and the variables x, y, x, and y are treated as independent. By collecting 
terms of equal powers of e, one obtains the leading solution involving (x, y), i.e., 
the analytical solution near each vortical spot. For the kth spot centered at 
(X k (t), Y k (t) ), the solution is (refs. 4 and 6) 



e -2-(0) 


r k 

2 6 

, 6 k (t) 


where 

6 k (t) = 6 k (0) t(t + t k )/t k ]V2 


(25) 


7 



is the effective core size and t^ = 6^(0 )^4v is the initial age of the vortical 
spot. Here the superscript (0) denotes "the leading term of.” The use of the 
similarity solution is explained in reference 4. Note that the solution given by 
equation (25) has only a "short range” effect because of exponential decay. On the 
other hand, the corresponding velocity is given by 


-1 

e 


<V* * V') 


[ E_1 'ViJv* - &)]/\ 


with 


-1 

e 


<v k 



) 


and 


(26a) 


(26b) 


\ = x - (x k /e) 


y k = y - <V«’ 


Here, e ^ (Vg)^ denotes the circumferential velocity around the kth vortical spot 
and does not decay exponentially but has a "long range" effect, i.e. , 


(v ) » 

1 e x 


2 nr. 


(27) 


This effect provides the coupling with the background flow. From an asymptotic 
analysis (refs. 4 and 6), the velocity of the center of the kth spot (X^, Y^) is 
established to be equal to the local velocity without the kth vortical spot, i.e., at 
(X k , Y k ), 


V 


V 


vv 


1 + 


U( W 


i + 


V< W 


A 'V ^ _<J _ A ^ A 

j + 2-j t (u.i + v j ) 


(28) 


j^k 


This result is in agreement with classical inviscid theory. The motion of the 
vortical spots^is, in turn, coupled with the variation of the background vorticity 
distribution £ that is 0(1) and is a function of the normal spatial variables. In 
subsection 3.2, the method for the solution of £(x,y,t) will be described. 
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3.2 Background Flow Field 

To derive the governing equations for the background flow field in the normal 
length scale, the details of the flow structure need to be filtered out in the small 
length scale eL. This procedure is accomplished by averaging the basic equations 
(eqs. (8) and (9)) over an area of, say, a square of side 2H where H is much 
larger than the core size eL but much smaller than L, i.e.. 


L >> H >> eL 


(29a) 


In the following equation, <f> is used to denote the average of f(x,y,x,y), i.e*/ 


2 A y+H)/e /*(x+H)/e 


<f> = 


Jfv-H) /e 


f (x,y,x,y) dx dy 


(29b) 


4H • / (y-H) /e •'(x-H) /e 


With f bounded for all (x, y) , the result is 


J&fr 

V3x> 


= 0 


df\ 

y9y/ 


= o 


(30a) 


For a function g(x,y) having no microstructure, i.e., independent of x and y, 
the result is 


<g> « g(x,y) 


(30b) 


Using equations (29b) and (30), the average of equations (8) and (9) yields the 
leading equations for the variation of background vorticity £(x,y,t) and the 
corresponding stream function ^(x,y,t); i.e.. 


& + 
d t dx 


{ [ u 0 (y) + u(x,y,t) + <e -1 u>]c} + |^-{ [v(x f y, t) + <e"'v>]?} 


- 1 -. 


+ 


[v + <e 


- 1 - , 
v>] 



0 


(31) 


and 



(32) 
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The initial and boundary conditions are, respectively 


C = o 


(t = 0) (33) 


and 


A/ 

= 0 (y = 0) (34a) 

8x 


VtJj * 0 


(|x| * y + 00 ) (34b) 


The averages <e~ u> and <e v> in equation (31) represent the bulk contributions 
of the velocity field induced by the moving vortical spots. It can be shown that 
these averages remain 0(1) with respect to e. (See the appendix.) It then follows 
that the grid size AL and the time step At can be selected independent of e or 
the core sizes 6^. 

The vorticity deviation £ is updated by the finite-difference equation 
(eq. (31)) using the two-step Lax-Wendrof f procedure (ref. 8). The fast Poisson 
solver (ref. 9) is then used to determine the corresponding stream function ij> from 
equation (32) for a finite computational domain. The computational domain and the 
approximate boundary conditions will be described in subsection 3.3. The trajec-^ 
tories of the vortical spots are then defined simultaneously by the integration of 
equation (28). 

Since the average has to be used only when the core size is much smaller than 
the grid size, the numerical method shall be discussed for the case of 6^ << AL. 

The integration domain size H is chosen to be of the order of the grid size (e.g. , 
H ■ AL/4). Therefore, H is much less than the normal length scale but is indepen- 
dent of the small core size. Whenever a grid point (x^, y j ) is far from a vortical 
spot centered at (X^, Y^) (e.g., if r^ » H, then r k > CH), the result is 


. -1- . - V 

<e V ~~2 

2irr, 


1 + 0 


H 


(35a! 


-1- r k (x i X k > 

<S V ' —2 

2irr, 


1 + O 


H 


(35b 


where 


= [(x. - X ) + (y. - Y ) ] 

c L i k i k J 


2 , 1/2 
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The difference between the average and the classical theory is actually less than 
0.6 percent for C = 2.5. Therefore, by adjusting C, the transition from the 
classical solution to the average can be made smaller than the error of the finite- 
difference method. Only in the region where r^ < CH do the averages <e u> and 
<e v> need to be evaluated (which are defined by two elementary line integrals, as 
noted in the appendix). 


3.3 Boundary Conditions on a Finite Computational Domain 

For the numerical solution of equations (31) to (34) for variations in the back- 
ground shear flow, one has to work with a finite computational domain D, i.e., 

|x| < L 1 and 0 < y < L 2< (See fig. 2.) Note that U Q (y) is a solution of the 
Navier-Stokes equation with a no-slip boundary condition. The boundary layer induced 
by the disturbed flow shall be ignored; hence, the vortical spot should be "far away" 
from the ground relative to the core size. Also, the use of the asymptotic solutions 
(eqs. (25) and (26)) for each core structure requires that the vortical spots be "far 
apart" from each other. Therefore, the present solution will be applicable only when 
each vortical spot is several core sizes above the ground and away from the other 
spots. In order to average over only one vortical spot, the requirement is that the 
distances from one spot to the ground and to another spot remain greater than twice 
the grid size, i.e., 2 AL. 

IV 

^ Note that the boundary condition ^ = 0 at y = 0 is fulfilled when the image 
5(x,-y) = -£(x,y) is introduced. The far-field boundary condition (eq. (34b)) has 
to be replaced by. appropriate conditions on the outer boundary of the domain D, 
i.e., Ixl = L, and y = L 0 . With p » (x^ + y2)1./2 denoting the distance from the 
origin, it is seen that when the disturbance velocity behaves as p , the behaviors 

<v /v iv 1 

of the perturbed stream function and vorticity £ are ii> * p and 

o* — n— 1 

C “ P , respectively. Since the flow field induced by the vortical spots is sym- 
metric with respect to the ground (y = 0), the variation ? of the background 
vorticity has to be antisymmetric with respect to y = 0 (with £ = 0 at 

y = 0) and approaches 0 near the outer boundary. Note that n = 2 and 3 for a 
single spot and for an equal and opposite pair of spots, respectively. Consequently, 

(V 

the boundary condition C = 0 can be imposed on the outer boundary (|x| = L 1 ) and 
on y = L 2 , with an error smaller than 0(.L^j for a single spot and smaller than 
0(L^ ) for a vortical pair of spots. On the boundary of D, the condition to be 
imposed is 


C - 0 


(36) 


Since the solution can be expressed in terms of a Poisson integral of £, the 
boundary data for ijj can then be evaluated by the far-field expansion of the Poisson 
integral (ref. 5). 

The leading two terms involve only the y- and xy-moments of i.e., at the 
boundary where y - L 2 or | x | = L 1 , the expansion is 


1 1 



( 37 ) 



2 2 1/2 -3 

where r = p = (x + y ) . The error of this approximation (eq. (37)) is 0{L. ). 

Since the error introduced by the finite-difference approximations for equations (31) 
and (32) are 0((AL) ), AL can be related to so that the errors are of the 

same order in the approximate boundary conditions and in the finite-difference ap- 
proximations. Using the boundary data from equation (37) and the condition in equa- 
tion (34a), equation (31) is solved by a fast Poisson solver (ref. 9) to obtain 
in D, and then the derivatives of t); yield u and v. 


4 RESULTS AND DISCUSSION 

A computer code based on the previous analysis was developed. Numerical results 
of the trajectories of vortical spots in a nonuniform shear flow are obtained. For 
the results reported here, the conditions = 8L, L 2 = 8L, and Ax = Ay = L/5 
have been chosen. The time step fulfills the condition that the Courant-Friedrichs- 
Lewy (CFL) number is less than 0.5. The velocity in the CFL number is the maximum of 
the resultant velocity in equation (31) for all grid points. Test cases were per- 
formed for larger L-j and b 2 and for smaller Ax and CFL number. In all these 
test cases, the changes in the results are insignificant. To gain a qualitative 
understanding of the interaction of vortical spots with a background shear flow, the 
case of a single vortical spot is studied initially. Then, the cases of a vortical 
pair simulating the trailing vortices are studied. 

Figure 3 shows the trajectories of a single concentrated decaying vortical spot 
of various strengths submerged in a background shear flow. The initial vertical 
position of the spot is at y = 1 , and the initial background shear flow is chosen as 
U Q (y) = 1 - e~ y . The results show that the vortical spots with positive circulation 
drift spanwise (x-direction) and upward (y-direction) , and the vortical spots with 
negative circulation drift spanwise and downward but eventually turn backward. This 
phenomenon is more pronounced as the strength of the vortical spot increases. 

To explain this phenomenon, the case of a single vortical spot with F > 0 is 
considered. The disturbed flow moves downward behind the spot (x < X) and upward 
ahead of the spot (x > X). For an initial background vorticity o)q with O)^(y) > 0 
(see eqs. (2) and (3)), the disturbed flow increases the vorticity behind the spot 
and decreases the vorticity ahead of it? i.e. , C > 0 for x < X and £ < 0 for 
x > X. The background vorticity variation £ in turn induces an upward motion of 
the vortical spot for F > 0. From similar arguments, it can be expected that the 
background vorticity variation £ will induce a downward motion of the vortical spot 
with F < 0. The reason that a vortical spot of negative strength turns around and 
drifts upstream as it gets closer to the ground can be attributed to the decrease of 
the contribution of the background shear flow to the forward velocity of the spot and 
to the increase of the induced velocity by the image of the vortical spot with 
respect to the ground (y - 0). It should be pointed out here once more that the 
vortical spot will drift horizontally when the background shear flow is either a 
uniform flow (o)q = 0) or a constant shear flow (idq = Constant), and there will be no 
change in the background flow (£ = 0). 
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The trailing- vortex wakes far downstream of an aircraft shall be simulated by 
a simple vortex pair whose vorticity distributions are concentrated and are centered 
at (iX^fO), Y^fO)) with strength ±T and effective vortical size 6(0). The goal 
of the following numerical examples is to simulate the interaction of the decaying 
trailing vortical pairs subjected to a crosswind (s panwise) ground shear. The back- 
ground shear flow used in the examples is again an exponential profile (eq. (2b)) 
unless otherwise specified. 

Figure 4 shows the relative trajectories, i.e., the variation of the vertical 
positions versus the horizontal spread of a vortex pair whose strengths are T = ±1 . 
Solid curves show the relative trajectories of a vortex pair descending without a 
background shear flow. Dashed-line curves show the corresponding trajectories of 
vortex pairs with an initial location at (±0.5, 3) and (±0.5, 4) under the influence 
of ground shear wind. The vertical position of the left vortex shows that it drifts 
to the left and continuously descends in time, whereas the right vortex descends 
initially until it reaches a minimum height and then drifts upward. Because x = 0 
has been set at the mean position of the vortical pair, this curve does not show the 
real spanwise position and does not see the eventual backward drift of the left 
vortical spot. This backward drift occurs when the background vorticity variation 
finally overcomes the downward velocity induced on the right vortex center by the 
left vortex. This effect is in contrast to the case when the background flow is 
uniform or of constant shear flow, for which the vortex spots will keep on drifting 
apart and downward. 

In order to find out when the shear- layer solution has to be used, the trajec- 
tories of a pair of vortices in a shear layer for different initial vortex heights 
are studied. The minimum height reached by the right vortical spot Y R m ^ n is 
plotted against the initial height of the vortical pair in figure 5. The height 
Yr min approaches an asymptotic value of 2.6 when Yj c (0) is greater than 7.0. This 
event means that when the vortical spots are above y = 7, they are far above the 
shear layer and the interaction with the shear layer is negligible. The correspond- 
ing trajectories (in real spanwise positions) of the vortical pair, starting at 
different heights Y^(-0) = 1, 2, 3, 4, and 5 in the shear layer, are displayed in 
figure 6 to show that the trajectories of the vortical spots are sensitive to the 
starting height (i.e., the altitude of the airplane relative to the thickness of the 
shear layer). 

The initial strength of the vortical spots also has a significant effect on the 
trajectories of decaying trailing vortices in a ground shear, as illustrated in fig- 
ure 7. The initial positions of the vortical spots are (X^(0), Yj c (0)) = (±0.5, 3.0). 
The vortical pairs of stronger strength (T = ±4) descend faster and drift less in 
horizontal (spanwise) distance as compared with the vortical pair of strength 
T = ±1, since the initial downward drift is dominated by their strength T. The 
vortical spot on the right side, with vorticity in an opposite sense to that of the 
shear flow, descends to a lower minimum elevation for larger T. 

Finally, the interaction effects of different background flows are studied. A 
comparison is made for a background exponential profile with a Blasius profile of the 
same thickness; i.e., the length of the profile is equal to L u ^/v. The initial 
positions of the vortex pair are (X^(0), Y^ ( 0 ) ) - (±0.5, 4.0) and their strengths 
are T = ±1 . Figure 8 demonstrates that the effect of the crosswind shear profile on 
the behavior of trajectories for the left vortex is small but the right vortex 
descends much less; i.e., it descends to a larger Y R m ^ n for the Blasius profile as 
compared with the exponential profile. 
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5 CONCLUDING REMARKS 

A computer code has been developed for the study of rotational flow fields 
having concentrations of large “magnitude vorticity in spots (small areas). The 
analysis takes into account the small scale of the flow field near a spot and the 
normal scale for the background flow. Using this two- length-scale analysis, the 
decay of the vortical spots (in small scale) is obtained analytically and the motion 
of the spot is coupled with transport of the background vorticity, which is analyzed 
by numerical solution of the Euler equation. The particular feature of the numerical 
method is that the grid size for the numerical analysis can be chosen independent of 
the size of the vortical core, and hence the time step is also independent of the 
core size. 

Numerical examples for a single vortical spot submerged in a background shear 
flow shows that the variation of the background shear flow induces an upward motion 
for the vortical spot with positive circulation and induces a downward motion for the 
vortical spot with negative circulation. The latter eventually reverses its direc- 
tion of horizontal drift when it gets close to the ground where the ground effect 
overcomes the smaller forward shear flow velocity. This interaction between the 
vortical spots and the background flow is absent if the background flow is of con- 
stant vorticity or is a uniform flow, for which the vortical spot will only drift 
horizontally. 

For a pair of vortical spots simulating trailing vortices in a cross wind, the 
result shows that both vortices will initially drift forward together, spread apart, 
and move downward. Gradually, under the effects of interaction, the redistribution 
of the background vorticity changes the trajectories. The vortex that is in the 
opposite sense to that of the background flow will reach a minimum height, reverse 
its downward drift, and turn upward, whereas the other vortex will eventually reverse 
its forward motion and turn backward. When the vortical pair is submerged in dif- 
ferent background flows, namely the Blasius and exponential profiles, the numerical 
examples show that the trajectories of the vortical pair are similar in nature but 
quantitatively different. 


NASA Langley Research Center 
Hampton, Virginia 23665-5225 
May 22, 1986 
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APPENDIX 


AVERAGE OP THE VELOCITY INDUCED BY kth VORTICAL SPOT <e _1 (u^i + v^j ) > 

It is necessary to derive the formula for the average velocity at a grid point 
(x^, y j ) that is close to the kth spot centered at (X k , Y^), with jx^ - X k | < H 

and Jyj - Y k J < H. (See fig. A1 . ) 

Because of the symmetry of the velocity field with respect to the vortical 
center (X, Y), the formula for <e u> shall be derived. The formula for <e v> 

is equal to that for <e” u> with x^ and X interchanged with y.. and Y, 

respectively. The restriction to be imposed is that 

H > X - x i > 0 H > Y - y^ > 0 (A1 ) 


Since u is asymmetric with respect to y = Y, the result is 

/. H y.2( Y-v . )-H . „ „ 

where 

r '= [( Xi + x' - X) 2 + (y + y' - Y) 2 ] 1/2 


-1- -r 

<e u> = — -z 


L dx ' J 


dy ' 


Yj + y -x 


/-« 


The integration can be carried out with respect to y ' and then a definite integral 
is obtained, i.e. , 



(A2) 


(A3) 
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where 


i 2 1 /$ 2 


and (R) is the exponential integral of R (ref. 10). 

From equation (A2) it can be shown that the integral with respect to R is 
finite; hence, <e -1 u> is 0(1) with respect to e. This result is also true for 
<e v>. 

Note that when R is greater than 2, ® 1 (R) is negligible and the integral in 

equation (A3) can be evaluated explicitly. Therefore, numerical evaluation of the 
integral in equation (A3) is needed only for r' < 26. 


(x. + x 1 , y. + y') 


v 


x., 

i’ 'J 


Figure A1 .- "Average" square domain 
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Figure 1 Physical model for motion and decay of vortical spots 
submerged in ground wind shear. 


y 


(-Lj, L 2 ) 

tL r 

4' 



D 



A?=-? 




j- Ground 

0, ?=0 



o 


I Image £(x, -y) = -£(x,y) ! 

i ——————— -i 


Figure 2.- Finite computational domain and boundary conditions 



r 

3 

2 



Figure 3.- Effect of strength on trajectoi 
vortical spot in shear flow. 
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Figure 7.- Effect of vortical-pair strength on trajectories 
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